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We provide a basis for the laser intensity dependence of the momentum distributions of electrons 
and ions arising from strong-field non-sequential double ionization (NSDI) at intensities in the range 
I = 1 - 6.5 x 10 14 W/an 2 . To do this we use a completely classical method introduced previously 
l|. Our calculated results reproduce the features of experimental observations at different laser 
intensities and depend on just two distinct categories of electon trajectories. 

PACS numbers: 32.80.Rm, 32.60.+i 



I. INTRODUCTION 

When atoms with two or more electrons are exposed 
to an ultra-short (10-600 femtoseconds) and ultra-strong 
(1 — 100 x 10 14 W/cm 2 ) laser pulse, two electrons can be 
ripped off from the nucleus. This strong-field double ion- 
ization process has been observed in all inert gas atoms 
IS El- tl and some molecules || . The remarkable aspect 
of this process is its strong electron pair correlation. Un- 
correlated electrons are ionized sequentially, each event 
well explained by single- ionization theory 6]. This se- 
quential ionization theory cannot explain a characteris- 
tic knee-like structure observed in the double-ionization 
yield (see Fig. [IJ as a function of laser intensity, so this 
process is usually known as nonsequential double ioniza- 
tion (NSDI). 

In addition to ion yield, experiments in the intensity 
regime of the knee-like structure have examined recoil 
ion momentum distributions in which the shape depends 
strongly on laser intensity 0, IS IS 13 13 13 El • Distri- 
butions for different atomic species (helium, neon and ar- 
gon) show different shape variations as a function of laser 
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FIG. 2: Experimental ion momentum distributions along the 
longitudinal direction for He 2+ |13| and Ar 2+ [t| at different 
laser intensities. As the laser intensity increases, the width 
of the distribution become broader and may develop into a 
double peak. We will comment on the experimental ion mo- 
mentum distributions of Ne 2+ at the end of this paper. 
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FIG. 1: Points trace out the knee-like structure of doubly- 
ionized ion yield in strong-field double ionization. The dashed 
curve (red) is inserted by hand to indicate the prediction made 
by sequential quantum tunneling theory 0. 
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intensity, as shown in Fig. |2 at low laser intensities there 
are more low-momentum or zero-momentum ions, giving 
a narrow distribution peaked at zero momentum, whereas 
at higher laser intensity the high- momentum ions become 
more dominant and the distribution becomes broader. 

A variety of theoretical initiatives have been pursued 
El 13 El IH 13 El with the goal of understand- 
ing NSDI effects including ion yield and final-state mo- 
menta. With the exception of the calculations made by 
the Taylor group that are applicable only to helium [15| , 
all have ad hoc elements that are not easy to justify 
fully, and none have successfully dealt with the intensity- 
dependence of the shapes of the experimental momentum 
distributions shown in Fig. PJJ 

Recently we pursued a fully classical theory 0] as an 
explanation for the electron correlation in NSDI and sug- 
gested that the major features of NSDI arise naturally 
without needing to engage quantum processes at all, not 
even tunnelling as the initial step in NSDI. In such a 
treatment the NSDI process is a consequence of both elec- 
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FIG. 3: Plots of electron energy versus time measured in laser cycles, showing four distinct stages of NSDI. We show a variety of 
histories of two-electron trajectory that lead to double ionization. One can identify four universal stages, which are embedded 
in these different histories of evolution. A trajectory is labelled as NZ and Z if both electrons moving away from the nucleus 
in the same and opposite directions respectively. See Q for detailed decription of the four stages. 



trons actively responding to the laser field and to each 
other and to the nuclear attraction, undergoing ordinary 
(but complicated) Hamiltonian dynamics throughout the 
duration of the laser pulse. As one consequence, double 
ionization can be easily seen to occur in four universal 
stages (see Fig. |3J). 

In this brief note we provide calculational evidence 
supporting a classical interpretation of the intensity- 
dependent changes in the shape of ion recoil momentum 
distributions that are shown in Fig. [3 The calculations 
suggest the active presence of two distinct categories of 
two-electron ejection from the atom, loosely speaking as- 
sociated with back to back ejection on the one hand and 
side by side ejection on the other hand. In the first case 
one can expect that little net momentum will be carried 
by the pair and in the second case the net electron mo- 
mentum can be substantial. For ease of description we 
have labelled these as the Z and NZ categories 0, , 
meaning "zero" and "non-zero" momentum to impart to 
the ion. 

Our approach will be to display the results of entirely 
classical calculations of the NSDI electron trajectories as 
evidence that the Z-NZ distinction exists, and then show 
that it explains the origin of the laser intensity depen- 
dence of the shape of the momentum distributions shown 
in Fig. H 



response of each electron is restricted to the axis of the 
laser polarization, will be adequate. All of the relevant 
experimental data is for linearly polarized laser light, 
which is the vast majority of available data. We also 
use the so-called Rochest er screened Coulomb potential 
[22| V(x) = —1/y/x 2 + 1 to model the nuclear binding 
potential and electron-electron repulsion potential. We 
have previously shown that the character of the NSDI 
process (recall Fig. |2Jl obtained using this Id model does 
not differ much from that of 3d In addition, Haan, 
et al. have shown that there are striking similarities in 
the dynamical evolution of the classical ensemble calcula- 
tions and the corresponding quantum 2-e wave functions 

mm 

Each electron pair is given a different set of initial con- 
ditions but each pair is given an initial energy of -2.24 
a.u., which corresponds to the energy of the two-electron 
quantum ground state. We use E(t) — £of(t) smut to 
denote the optical field, where So is the maximum field 
strength, f(t) is the envelope function. We take an 8- 
cycle (~ 25 fs) sinusoidal laser pulse with the wavelength 
780 nm (frequency u = 0.0584 a.u.) and a trapezoidal 
envelope with 2-cycle turn-on and turn-off. 

The response of the two electrons can be obtained by 
solving two coupled Newtonian equations of motion: 



II. CLASSICAL ENSEMBLE METHOD 

We use the classical ensemble method that has been 
described in detail by Panfili, et al. |21j. The idea is to 
use a large number of classical electron pairs (100,000 - 
500,000 members) to model the dynamical evolution of a 
two-electron system exposed to an intense laser field. 

We are interested only in the forward-backward as- 
pects of the recoil ion momentum, what is measured as 
the longitudinal momentum (i.e., along the axis of the 
laser polarization) . This means that the one-dimensional 
aligned-electron approximation ( AEA) , in which the 
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Here X\ and x-i are the positions of two electrons, and 
clearly the evolution of each electron pair can be visual- 
ized as a distinct trajectory in the x\ — x 2 plane. In 
this treatment all Coulombic forces and laser-electron 
forces are acting continuously and simultaneously, and 
both electrons are continually dynamically active. 
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FIG. 4: End of pulse momentum distributions of doubly ion- 
ized electrons at various intensities. These are density plots 
in log scale, in which regions with darker color have relatively 
higher densities of two-electron trajectories. The box size in 
the figures is much smaller than the computational zone and 
was chosen for ease of viewing. It does not represent a bound- 
ary for the electrons. Notice that there are no trajectories 
close to the two axes due to our definition of double ioniza- 
tion, in which an electron is considered ionized if its kinetic 
energy is greater than its nuclear binding potential energy 
. If the definition is modified, the fine details of these dis- 
tributions will change but the main features will remain the 
same. 



III. CALCULATED MOMENTUM 
DISTRIBUTIONS OF IONS AND ELECTRONS 

Calculations based on the ensemble method provide a 
wealth of "experimental" information at every intensity. 
Fig. 0] shows the two-electron momentum plane at the 
end of the laser pulse for each of six different laser inten- 
sities. The planes contains dense distributions of doubly 
ionized electron-pair points. Electron pairs with oppo- 
site momenta (Z category) produce points in the second 
and fourth quadrants, and electron pairs moving together 
(NZ category) provide the points in the first and third 
quadrants. These points yield theoretical ion momentum 



distributions for the Z and NZ pairs, and these distribu- 
tions are found to have different shapes, and also differ- 
ent behaviors as a function of intensity. The calculated 
shapes are shown with their intensity dependences in Fig. 
[5] In both categories the shapes broaden with increasing 
intensity, more so in the NZ case. More significantly, in 
both categories the peak of the distribution drops with 
increasing intensity, but this effect is much stronger for 
the Z category. Analysis of the time evolution of Z and 
NZ trajectories, and their relative ability to produce the 
NSDI "jets" that have been discussed previously p^l23 |. 
will be reported separately. 

The combined result of broadening and peak-height 
decrease is a predicted substantial shift in probability in 
favor of the NZ category. As shown in Fig. H3 the rel- 
ative probability P(Z) - P(NZ) is predicted to change 
sign as the intensity exceeds a few times 10 14 watts/cm 2 , 
first strongly favoring the Z category at low intensity and 
then above about J~5x 10 watts/cm 2 strongly favor- 
ing the NZ category. This correlates with the change in 
appearance of the panels in Fig. 01 and we can see that 
this Z — > NZ shift is adequate to explain the shape varia- 
tion with increasing intensity in the experimental curves 
shown in Fig. [21 Additional quantitative comparisons 
can be obtained if the Z and NZ data are combined to 
make calculated total-momentum distributions. 

Fig.[7]plots the calculated ion momentum distributions 
obtained by combining Z and NZ contributions for inten- 
sities in the range of / = 1 — 6.5 x 10 14 W/cm 2 , similar to 
the range in some of the experiments. Details of the dis- 
tributions depend to a degree on the definition of double 
ionization, as explained already in Fig. ^ but the main 
characteristics are quite stable. By comparing our calcu- 
lated results with those measured experimentally (recall 
the data in Fig. |2J), we find the following similarities in 
the laser intensity dependence of the shape of the mo- 
mentum distribution: 

(1) The width of the distribution increases with the 
laser intensity. The helium experiment shows that the 
base width of the distribution changes from 8 a.u. at 
3.5 x 10 14 W/cm 2 to 12 a.u. at 12.5 x 10 14 W/cm 2 , and 
the argon experiment shows that the base width of the 
distribution changes from 10 a.u. at 2.0 x 10 14 W/cm 2 
to more than 12 a.u. at 3.75 x 10 14 VF/cm 2 . Sim- 
ilarly, our calculated distribution has the base width 
of 6 a.u. at / = 1.0 x 10 14 I^/cm 2 . This width is 
raised more than 50% trajectories to about 10 a.u. at 
/ = 6.5 x 10 14 VF/ cm 2 . 

(2) More low or zero momentum ions are observed at low 
laser intensities, whereas more high or nonzero momen- 
tum ions are observed at high laser intensities. Such a 
transition is clearly illustrated both in the experiment 
and our calculation, in which the side peaks become 
higher and wider as the laser intensity increases. 

One can see from the calculated two-electron momen- 
tum distributions that four NZ regions near one axis but 
relatively far from the origin have the highest density of 
NZ trajectories, and two Z regions near the origin have 
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FIG. 5: The recoil ion momentum distribution of Z and NZ trajectories at different laser intensities. Each distribution is 
calculated by first projecting the two-electron momentum points in Fig. [I] onto the diagonal axis. The raw data projections 
are then grouped in bins with width 0.3 a.u., and their number is divided by the total number of NSDI points to give the 
probabilities plotted here. 
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FIG. 6: The fractional yield of Z and NZ trajectories at dif- 
ferent laser intensities. Note that P(Z) and P(NZ) indicate 
the probability of yielding Z and NZ trajectories respectively. 



the highest density of Z trajectories. The presence of 
these densely populated regions indicates that the mo- 
menta of two electrons are correlated. We find that our 
classical calculations show similarities with the quantum 
results of Lein, et al. j2||, obtained by solving the time- 
dependent Schrodinger equation for the same model with 
similar parameters. The locations of the highest density 
two-electron trajectories are in good agreement with the 
quantum calculations. This is a further confirmation that 
the dynamics of two electrons in a strong laser field has 
definite classical attributes. 



IV. SUMMARY 

We have used a completely classical method to inves- 
tigate the laser intensity dependence of the momentum 
distributions of electrons and ions. Based on the final 



5 




Ion Momentum (a.u.) 



FIG. 7: Momentum distribution of the recoil ion at various 
intensities. To reveal the differences at different laser inten- 
sity, the calculated data points, which are grouped in bins 
with width 0.3 a.u., are connected with line. 

momenta of two electrons, we have identified two dis- 
tinct categories of electron trajectories, Z and NZ. We 
found that the shapes of the ion momentum distribu- 
tions of the Z and NZ trajectories are very different, in 
which the Z and NZ trajectories produce single-peak and 
double-peak structures respectively. Furthermore, the 



height and width of their distributions depend on the 
laser intensity. Thus, these two types of trajectories act 
as a basis to give different shapes of the ion momentum 
distributions at different laser intensities. Our calculated 
results showed that the Z trajectories are more favorable 
at low laser intensities, whereas the NZ trajectories are 
more favorable at high laser intensities. This Z — > NZ 
shift has been demonstrated in the helium and argon ex- 
periments. 

Currently, the available experimental ion momentum 
distributions of neon show a double-peak structure, and 
its width decreases with decreasing laser intensity within 
the NSDI regime (see 0,0)0)- Since this shape and 
its variation are consistent with the behavior of the NZ 
distribution, we suggest that these neon data are in the 
intensity regime that favors the production of NZ trajec- 
tories 26]. By lowering the laser intensity further, our 
analysis suggests that experiment on neon could observe 
a sngle-peak momentum distribution that has dominant 
Z ions. 
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